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I. INTRODUCTION 

It is well known that the Euler-Heisenberg effective 
Lagrangian in quantum electrodynamics (QED) is pre- 
sented by an asymptotic series jli]. Since the asymp- 
totic series is divergent its appHcations are limited to a 
weak field approximation. For the problems related to a 
strong field background one has to apply the complete ex- 
pressions for the effective action. The effective action at 
one- and two-loop level for special electromagnetic back- 
grounds have been studied in 0, Q (see as well and 
refs. therein). Recently the analytic series representation 
for the one-loop effective action of QED has been con- 
structed [3 from Schwinger's integral expression for the 
effective action The knowledge of an explicit analyti- 
cal expression of the effective action is useful for studying 
quantum effects like photon splitting and pho ton prop- 
agation in electric and magnetic fields l9l.llCt [lll | . An- 
other motivation comes from the theoretically predicted 
existence of strong magnetic fields of magnitude B above 

2 

the critical value Bcr = — 5 • 10 G in the cores of 

e 

some neutron stars 12]. 

In the present paper we consider light propagation ef- 
fects in the constant magnetic and electric field in the 
framework of the effective action approach proposed in 
jll'l applying the series representation for the one-loop 
QED effective action |5|. In Section^ we present some 
general formulae relevant to the effective action approach 
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In sections IlIII and ITVI we consider photon propaga- 
tion in the pure magnetic and pure electric external field 
correspondingly. Birefringence effect in the magnetic and 
electric field is considered in sections lVllVIl Some techni- 
cal details of calculation and useful formulae are collected 
in the appendix. 



II. LIGHT CONE CONDITION IN THE 
EFFECTIVE ACTION FORMALISM 

The integral representation for the one-loop effective 
action of QED for slowly varying background fields has 
been known for a long time [Sj 

A/: = -^ah / - coth(eai) cot(e6i)e^'" *, (1) 

where a, h are gauge invariant variables corresponding to 
magnetic and electric fields in a Lorentz frame where the 
both field vectors are parallel 

F'^" = \^^'''^Ff,„. (2) 

We start from the series representation for the one- 
loop effective Lagrangian which is more suitable in 
physical problems when we need explicit analytical ex- 
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pression of the effective Lagrangian 
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where Cci is the classical QED Lagrangian. 

We will follow the effective action approach 11] to 
study the light propagation effects in non-trivial vacua. 
The light cone condition for a soft photon propagating in 
the homogeneous electromagnetic background (neglect- 
ing the back-reaction effects) can be derived from the 
equation of motion corresponding to the full effective ac- 
tion. The final expression for the light cone condition is 
the following 







(«,k), 
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where is a unit polarization vector of the photon, and 
is defined by (in the Lorcnz gauge fc^e^ = 0) 
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The matrix -/V/^^ provides the quantum corrections to the 
classical light cone condition fc^ = 0. The light cone con- 
dition can be also transformed into another form which 
is more suitable for calculation of the light velocity aver- 
aged over polarizations. After averaging over the polar- 
ization states of the photon the light cone condition can 
be rewritten as |0 



= 1 - Q < T'"' > fc^fc,, 



(6) 



here, < T'"' > is a vacuum expectation value (VEV) of 
the electromagnetic energy-momentum tensor, and Q is 
an effective action charge defined in terms of the effec- 
tive Lagrangian and its derivatives jllj | . For convenience 
we give the definition of the effective action charge Q in 
terms of gauge invariant variables a, b 



Q = [dlC + dlC)[C{dlC + dlC) + 

-^L^iadaC - bdbCf - adaCdlC - bd^Cdlc) ~\{7) 
+ b"^ J 



As can be seen from the light cone condition, Eq. lO, 
the light velocity is determined by the product of the 
effective action charge Q and the VEV of the energy- 
momentum tensor < T^i, >. It is convenient to introduce 
the following notations for the combinations of special 
functions 



g{x) — ci(x) cos X -I- si(x) sin a;, 
h(x) = i [e^Ei(-a;) -I- e-^Ei(a;) 



(8) 



The properties of the special functions are given in 
H^ . Using Eq. (j2Jl for the one- loop correction (5£, one can 
write down explicit expressions for the derivative terms 
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a na 2 \ , /-\ 

— coth a H ; — csch a n[z) 

n b ' 



-h'{z) coth a 



disc = [—coth /3 + .g(z) coth /?(- 
ijT'^ ^-^ Vna \ 



47^3 z — i ^^fj^ 

n— 1 



-csch /9 



a 



-c&ch^ I3g' (z) 



27r 



coth a ]csch^ ahlz) 

b^ b 



„, - ^ e^a sr-^ r 1 , , , /' 

o^oC = ^— g 1^ — — cotha -I- /i(z) cotha(^- 



2&3 



— -- — csch a 
b-^ 

27r 2&7r2n 



2am: 
63 e 



csch^a/i'(z) 



coth (3 ) csch^ (3g{z) 



(9) 



where 
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The vacuum averaged value of the energy-momentum 
tensor can be found from the effective Lagrangian as in 
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here, T^j^ is the classical energy-momentum tensor 
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(12) 
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In the following we will apply the above general equations 
to special cases of light propagation in homogeneous pure 
magnetic and electric background. 



III. LIGHT VELOCITY IN MAGNETIC FIELD 



Light propagation in the presence of a constant mag- 
netic field in one-loop approximation was first studied in 
0,11,01 and later by many others (see refs. in Q). The 
light velocity in the weak magnetic field was well stud- 
ied using the Euler-Heisenberg effective Lagrangian 
(weak-field limit of the complete one-loop effective QED 
Lagrangian) 



2m f ea \^ 
C-ci H J- 1 



^ (271 2) (2n + 3) (2?i + 4) V 



(13) 



where m is the electron mass and _B„ is the Bernoulli 
number. Clearly the series HI3|I is an asymptotic series 
which is divergent, so that for ultra strong magnetic field 
problems one can not apply it for computational pur- 
poses. For instance, for the critical value Bcr of the mag- 
netic field, which supposed to be in the cores of some 
neutron stars, the series in the Euler-Heisenberg effective 
Lagrangian Ijf 3|l gives alternative in sign values when we 
include more than three terms in the truncated series. 

In this section we evaluate the light velocity using 
the series representation for the one-loop effective La- 
grangian, Eq. (O, in weak and strong magnetic field re- 
gions and derive asymptotic formulae in the ultra strong 
field regime. 

Without loss of generality we choose the magnetic 
field to be directed along the z-axis, so that B = (0, 0, a). 
The non- vanishing components of the classical and quan- 
tum energy-momentum tensor arc the following 
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< T^^ >-- 

= — SCrn - adaSC. 



(14) 



Hereafter in this section the functions daSC, dbSC, d^5C, 
&l5C are taken in the limit of vanishing electric field, 
6^0, 



dadC = -rG{a) 



dtS£ = 0, 
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The function G{a) defines the one-loop contribution to 
the effective Lagrangian with a pure magnetic back- 
ground and can be written in terms of the generalized 
Hurvitz ^-function. From Eq. one can find imme- 
diately that the velocity of light pr opag ating along the 
magnetic field remains unchanged |9lllf| whereas for the 
photon moving in the plane orthogonal to the magnetic 
field one has the phase velocity given by 
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The VEVs of the energy-momentum tensor can be writ- 
ten in terms of the main function G(a) 
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It should be noted that after taking the limit 6 — > in 
the general formulae © we have obtained a non-trivial 

contribution from the non-vanishing term d^C . This 

contribution provides the leading decrease factor in the 
dependency of the light velocity in the strong field region. 
After simplification the final expression for the light ve- 
locity takes the form 
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To compare our results with the results obtained in past 
we consider the light velocity in weak field approxima- 
tion which can be obtained using the truncated one-loop 
effective Euler-Heisenberg Lagrangian 0, ^| 
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(20) 



4 



0,2 



0.4 



0.6 



0,999998 
0,999996 
0,999994 
0,999992 
0,99999 
OS 




Notice that the trace anomaly < T''^ > of the energy- 
momentum plays a crucial role in providing the upper 

' - >. < 1 



bound for the factor Q{< T"" > 



2 < 

until a reaches the value of the Landau pole. Our for- 
mula improves on previously known expressions for the 
light velocity 0| and can be applied to ultra strong 
magnetic fields. 



IV. LIGHT VELOCITY IN HOMOGENEOUS 
ELECTRIC BACKGROUND 



FIG. 1: Light velocity in the plane orthogonal to the magnetic 
field: (i) exact one-loop approximated result; (ii) result in the 
weak field approximation. The magnetic field strength a is 
measured in units of electron mass m^. 



Let us consider an homogeneous electric field E = 
(0,0,6) directed along the z-axis. The non-vanishing 
components of the energy-momentum tensor and its 
VEVs are 



For the photon propagating in the plane orthogonal to 
the magnetic field one can derive as in the following 
expression for the velocity 
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1/137.036 is a fine structure con- 



stant. The corresponding light velocities obtained from 
the exact one-loop effective Lagrangian Q and one-loop 
Euler-Heisenberg Lagrangian (|20(l are plotted in Fig. 1. 

Now, let us consider light velocity in the strong field 
region. From the asymptotic behaviour of the main spe- 
cial functions ci(x), si(x) 'T3,'T4j one can find the asymp- 
totic formula for the function G{a) 
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With this one results in the following asymptotic 
equation for the light velocity in a strong magnetic field 
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The last term in the denominator in Eq. H23|) gives a lin- 
ear decrease of light velocity ( for strong magnetic fields 
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urn 



in qualitative agreement with results in 



As in the case of a pure magnetic background the veloc- 
ity of light propagating along the electric field remains 
unchanged and in the directions orthogonal to the elec- 
tric field we have for the light velocity the same equation 
(16). The formula for the light velocity derived from the 
one-loop Euler-Heisenberg effective Lagrangian coincides 
with one expressed by Eq. (19) with the replacement 
B ^ E. The explicit expressions for the effective action 
charge Q and VEVs of the components < r°" > , < > 
can be derived in a similar manner as in the previous 
section 
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The function H{b) is defined by 
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FIG. 2: Light velocity in the orthogonal plane to the electric 
field: (i) the exact one-loop approximated result; (ii) the re- 
sult obtained in weak field approximation. The electric field 
strength b is given in units of electron mass m^. 



In the asymptotic region approaching the Landau pole 
from the left the light velocity squared takes negative 
values. This shows that the Landau pole problem is still 
present. 



V. BIREFRINGENCE IN MAGNETIC FIELDS 

In this section we consider the birefringence effect 
with a complete one- loop effective Lagrangian Q. Let us 
choose the direction of the homogeneous magnetic field 

along the z-axis and the wave vector k in the a;Oz-plane 

— [v, shiO, 0, COS0), 
= (0, 0, 1, 0), 
e|[ = (0, -cos 6*, 0, sine), 



(30) 



After simplifications we find the following equation 
for light velocity in the orthogonal plane to the electric 
field 



2b + dbSC - bdlSC 
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where 9 is the angle between the vectors B, k and e'[, 
are polarization vectors in orthogonal and parallel direc- 
tions to B,k-plane correspondingly. 

To consider birefringence in the magnetic field we will 
start from the light cone condition Eq. Q . The non-zero 
components of the matrix Af^^ are 

Mil = a'dlC 



-"'^03 



^(9 C 

y 



(31) 



The exact values of the light velocity v in the weak field 
region b < 1 are plotted in Fig. 2 and can be compared 
with ones obtained from the one-loop Euler-Heisenberg 
Lagrangian if^ . 

Of course, for a strong electric field of magnitude 

above the critical value Ecr = pair creation, which 

e 

we have not taken into account, becomes important, so 
in the strong field region one can only consider the light 
velocity dependence on the electric field b formally. In 
that case the light velocity dependence is determined by 
the asymptotic behaviour of the function H(b) 



H{b) = _!^(ln4 + ci) 
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where we keep only leading terms. 
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For a strong electric field one has the following asymp- 
totic formula 
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With this the light cone condition can be written as 
follows 



(32) 



For the _L, || -modes the last equation is simplified to 

{vl-l)da,C = a^sin^edlC, (33) 
{vj-l)da,C = a^sin^ed^Cvf. (34) 

With this one can obtain the light velocities for _L, || po- 
larization modes 
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FIG. 3; The upper curve, (ii), corresponds to the ||-mode and 
the lower one, (i), corresponds to the ±-mode. This can be 
compared with the light velocity averaged over polarizations, 
curve (iii). 



For the special case ^ — 1^ the light velocities for (_L, ID- 
modes are plotted in Fig. 3. 

The estimated values for the velocities t)j_ y can be 
compared with ones obtained in Q. Tsai and Erber de- 
rived the explicit expressions for the refraction indices for 
_L, II modes [9 = 7r/2) 



where h 



^±,11 = 1 
m 



■ o{a 



(39) 



2eH 



(H coincides with our notation a) and 

77 functions are expressed in terms of special functions. 
This formula was derived approximately by neglecting 
terms of order higher than a. So, the Tsai-Erber' formula 
is not valid for large magnetic fields a >> 1. For weak 
fields our results are in a good quantative agreement with 
Tsai-Erber' results [tit-e) (see Table I). 
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0.01014 
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0.004725 
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0.00482 


0.004801 


0.00270 


0.002708 
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0.00309 


0.003082 


0.00175 


0.001748 



TABLE L Comparison of the Tsai and Erber results with 
our results for the functions ?7||,x. 

For the strong field region a >> 1 one can find the 
asymptotic formulae for the light velocities 
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In the asymptotic limit a 
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One can also find the corresponding refraction indices 
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For instance, keeping only the leading linear term in H4()|l 
one can approximate the refraction index n\\ by 



1 



T^fr ' (43) 

in agreement with other authors [lol ITsj l. Notice also 
that the birefringence effect was considered recently in 
where the effective action approach Till had been 
applied to a one-loop part of the Euler-Heisenberg effec- 
tive Lagrangian (|20|l . 



VI. BIREFRINGENCE IN ELECTRIC FIELDS 



Let us choose the direction of the electric field along 
the 2:-axis, so that 



i^'"'^ = except = |E| = 5, 
F^''' = except ^ ^ 



(44) 



The photon momentum k and polarization vector e are 
the same as ones considered for the pure magnetic case. 
The non-zero components of the matrix M^^f are given 

by 



MS! = -b^dlC 

f03 



(45) 



With this the light cone condition (0J can be rewritten 
as follows 



-2M°!vte''es. 
For _L, II -modes the last equation gives 

{vl-\)d^C = b^sm^edlCvl- 
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nally 
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useful identities 
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FIG. 4: Light velocity in the electric field for the ||-mode, (i), 
and for the ±-mode, (ii). 



where derivatives of the 5C are taken at the limit a = 0. 
For electric fields of magnitude below the critical value 
6 < 1 one has a small deviation of the light velocity for 
_L, 1 1 -modes from the velocity averaged over polarizations 
as it is pictured in Fig. 4. 



VII. CONCLUSION 



In this paper we have demonstrated the applications 
of the analytic series representation for the one-loop ef- 
fective action of QED proposed recently [5|. We have 
considered the light propagation effects in constant elec- 
tric and magnetic fields and demonstrated the advantage 
of our effective action approach in comparison with the 
approach using the asymptotic series of Euler-Heisenberg 
effective Lagrangian (weak-field approximation of the 
complete one-loop effective Lagrangian) in the case of 
strong and ultra-strong magnetic fields. The analytical 
asymptotic formulae for light speed and the refractive in- 
dices for arbitrarily strong magnetic fields have been ob- 
tained. Even though our analysis is restricted to one-loop 
consideration, we hope that our results describe qualita- 
tively the physical effects due to the very small value of 
the fine structure constant, and can provide good quanti- 
tative results in astrophysical problems concerning strong 
magnetic fields. The generalization to two-loop order is 
of further interest. 

APPENDIX 

A. Notations and useful equations 

Our conventions for the metric tensor (^^^ and anti- 
symmetric tensor e^'^'"^ are the following 



g0123 
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We use Greek letters for 4-dimensional space-time coor- 
dinates and Latin letters for space indices only. One has 
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Relations between the variables a,b and x,y, and corre- 
sponding derivatives are given by 
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d,y = —{abdl-abdl + {a^-b')dab)- 
\\{ia% - b'^^da + (3a62 - a^)dbl 



(52) 



Using the change of variables (x, y) — > (a, b) one can 
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express the matrix M^^ in terms of variables (a, b) 
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2abF^''F^p + {a^ - b^){F^-'F^p + F'^'' F^0)]dabC 
^[(3a62 - a3)F'^"^i^„^ - {3ab^ - a^)F'^'' F^p 

iiaH - b^){F^''F^t, + F^'''F^p)]daC 
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(Sab^ - a^){F'^-F^p + F^'' F^p)]dbC. 



(53) 
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For completeness we write down the expressions for 
the energy-momentum tensor TL^ and effective charge Q 
in terms of the cc, y variables [llj 



dS 



S = J d'^x^—gC, 
< T^"^ >= -T^"'^^C + g'^-'iC - xd^^C - ydyC), 

-ryd,y)L + \(dl + dl)L{\ - xd^ - vdy)L) '\ (54) 
B. Asymptotic formulae 

The asymptotic properties of standard special func- 
tions are given in 14]. Keeping the leading terms in 
series expansions we obtain the leading terms in asymp- 
totic expansion of the main special functions: 
in a; — > oo region 



e-="Ei(x) 
e^Ei(-x) 

g{x) 

h{x) 



1 

X 



(55) 



in a; ^ region 



Ei(x) 
Ei(-x) 
ci(x) 

si(x) 

g{x) 

h{x) 



7 -I- Inx, 
7 -I- Inx, 
7 -f Inx, 

TT 

7 4- Inx, 
7 -|- Inx. 



(56) 



The asymptotic formulae for the functions G{a), Eq. 
(|22|l . and H{b), Eq. (|28|) . can be derived in several ways, 
either from Schwinger's original integral formula for the 
one-loop effective Lagrangian using the generalized Hur- 
witz C-function, or using the asymptotic formula for g{x) 



and the series representation Eq. Consider the inte- 
gral representations for the main functions 



9{x) 
h{x) 



1 + ^2-^*' 
dt. 











te- 



(57) 



Substituting these expressions in Eqs. IjlGI I26() and per- 
forming the summation one finds the integral represen- 
tations 



G{a) = - 
H{b) = 



dtt 



■Li2(e-"^*), 



dtt -IHBlf. 

^^Li2(e e*. *) 



(58) 



from which the asymptotic formulae for the functions 
G(a), H{b) can be easily derived. The values of the func- 
tions G{a),H(b) obtained from the series representation 
and from the asymptotic formulae for large a, 5 ^ 1 can 
be viewed in Table II. 



a 


G series 


G asym 


b 


H series 


Hasym 


10^ 


-0.37815 


-0.316198 


10^ 


+ 0.122089 


+ 0.081690 


10^ 


-2.38586 


-2.385448 


10^ 


- 2.082565 


- 2.082954 


10^ 


-5.71865 


-5.718648 


10^ 


- 5.653113 


- 5.653116 


lO"* 


-9.42814 


-9.428142 


10'' 


- 9.41787 


- 9.417870 


10^ 


-13.2046 


-13.204571 


10^ 


-13.20317 


-13.203169 


10'' 


-16.99072 


-16.990717 


lO'^ 


-16.99054 


-16.9905396 



TABLE 11: Comparison of the results obtained from series 
representation with the results derived from the asymptotic 
formulae 
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